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The A + A — > 0, B + B — > process with exclusion be- 
tween the different kinds is investigated here numerically. Be- 
fore treating this model explicitly, we study the generalized 
Domany-Kinzel cellular automaton model of Hinrichsen on 
the line of the parameter space where only compact clusters 
can grow. The simplest version is treated with two absorbing 
phases in addition to the active one. The two kinds of kinks 
which arise in this case do not react, leading to kinetics differ- 
ing from standard annihilating random walk of two species. 
Time dependent simulations are presented here to illustrate 
the differences caused by exclusion in the scaling properties of 
usually discussed characteristic quantities. The dependence 
on the density and composition of the initial state is most 
apparent. Making use of the parallelism between this pro- 
cess and directed percolation limited by a reflecting parabolic 
surface we argue that the two kinds of kinks exert marginal 
perturbation on each other leading to deviations from stan- 
dard annihilating random walk behavior. 



I. INTRODUCTION 

Non-universal dynamical behaviour seems to be a con- 
troversial issue in non-equilibrium models. An outstand- 
ing example is the debated behaviour of systems exhibit- 
ing infinitely many absorbing states There is no 
analytic treatment up to now; argumentation of various 
authors, in most of the cases, is based on simulation re- 
sults. Despite intensive study, the critical behaviour of 
such systems is poorly understood, non-universality re- 
mains an unresolved problem and even scaling behaviour 
is questioned. Roughly speaking, in these coupled pro- 
cesses the 'primary' particles follow a branching and anni- 
hilating random walk while the other species just provide 
a slowly changing environment that effects the branching 
rates of the primaries. The spreading exponents of the 
primaries depend on the initial conditions of the environ- 
ment. 

A possible way which might lead to a deeper under- 
standing of the mechanism behind non-universal spread- 
ing could be the study of simpler coupled systems. Per- 
haps the simplest case is the coupled annihilating random 
walk of two species ( A + A ^ %, B + B ^ $). Naively 
one would expect that this could be described by the 
exactly solved field theory of the A + A — ► process ||] 



(ARW). In one dimension, however, the situation is more 
subtle than in higher dimensions. Particles of different 
type can block the motion of each other. The difference 
between one and two dimensions has been found to give 
rise to different phase diagrams in the case of the general 
epidemic model [Q. The question now arises how rele- 
vant the exclusion perturbation caused by this blocking 
mechanism is to a fixed point of the kind determined in 
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An other motivation of this study originates from the 
investigations of Hinrichsen j^], who found, by simula- 
tions, a str ang e scaling behaviour in some special case of 
his model |2l[] ,for which, however, an explanation is still 
lacking. In section |0| Hinrichsen's model will be intro- 
duced. It is easy to see, that the kinks in this model at the 
symmetry point corresponding to the compact directed 
percolation point of the Domany-Kinzel automaton, ex- 
hibit the process described above. In sections III and IV 



we present our high precision time dependent simulation 
results from random and seed initial conditions. In sec- 
tion |v| these results are compared with those obtained 
by rigid (i.e. parabolic) boundaries. We further investi- 



gate this an alogy on the mean-field level in section VI 



while section VII is devoted to results in the explicit two- 
species Annihilating Random Walk model with exclusion 
(ARW 2e). We summarise our numerical results in section 
VIII and give an outlook toward iV-species generalisation 
A qualitative description of the behaviour of Hin- 



IX 



richsen's model outside the symmetry point on the line 
of compactness is presented in section [x] and finally in 



section 



XI we summarise and discuss our results. 



II. THE GENERALISED DOMANY-KINZEL SCA 

The Domany-Kinzel (DK) stochastic cellular automa- 
ton (SCA) § is one of the simplest models which show a 
non-equilibrium phase transition into an absorbing state. 
This one dimensional SCA is defined on a ring with two 
states '1' and '0' with the following rule of update: 
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where at t + 1 the probability of 'l'-s is shown . 

In the plane of the parameters (p, q), the phase diagram 
of the DK SCA is as follows. A line of critical points 
separates the active phase ( with a finite concentration 
of 'l'-s) from the absorbing (vacuum) phase (with zero 
steady state density of 'l'-s). This continuous transition 
belongs to the universality class of directed percolation 
(DP) @. The endpoint of this line (q = 1, p = 1/2) 
describes a transition, however, outside the DP class; it 
corresponds to compact directed percolation. Here the 
model exhibits Ising symmetry and can be solved exactly 



and the phase transition line now belongs to the parity 
conserving (PC) universality class. This class has been 
studied by many authors as the first exception from the 
robust DP class 0,P,&0Jl|,|2fll . 
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FIG. 1. a) .The phase diagram of the DK SCA and b). the 
corresponding phase diagram in case of the simplest version 
of the GDK cellular automaton model by Hinrichsen 



In 1997 Hinrichsen Q introduced a generalized version 
of the DK model including more than one symmetric in- 
active states (71, 72, ...) and one active state (A). The 
motivation for this study was to look for a possible change 
in the universality class of the line separating the active 
and passive steady states. This generalized DK model 
(GDK in the following), in its simplest form with two 
absorbing states 71 and 12 has been defined by the rules 
given below: 
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The geometry of updating is the same as in the case of 
the DK SCA. It has been shown by simulation || that the 
phase diagram which emerges is similar to that in the DK 
SCA : an active phase is separated from an inactive one 
by a line of continuous phase transitions. The inactive 
phase, however is symmetrically degenerated (II or 12) 
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FIG. 2. Evolution from a random initial state in the gen- 
eralised Domany-Kinzel SCA (GDK) on the line of compact- 
ness, at q = 1, p = 0.5. Red: 71, blue: 72, white: A 

The phase diagram exhibited on Figl.b) shows that 
the line of PC transitions ends at q = 1, p = 1/2, a point 
which corresponds the Ising symmetry point of the DK 
automaton. The primary aim of the present work is to 
investigate the scaling properties of GDK at this point, 
which will be called CDP2 transition point. A typical 
time evolution of the GDK model at this special point 
when starting from a random initial arrangement of 71- 
s, 72-s and ^4-s is shown on Fig. || Here active islands 
can be spatially extended, thus three kinds of compact 
clusters can grow. Nevertheless only the 71 and 72 phases 
are ^-symmetric while the active phase plays a special 
role. (The situation is different from a 3-states Potts 
model with Glauber kinetics ). 

It is well known that the CDP process in Id is equiv- 
alent to an annihilating random walk process of kinks 
separating compact domains of 0-s and 1-s. In the 
model investigated here two types of kinks can be de- 
fined, namely kink Kl between domains A — 71 (and 
II — A) and kink K2 between neighbouring A — 72-s (and 
72 — As) . The rules of the model inhibit occurrence of 
kinks between domains of absorbing phases, i.e. between 
71-71 and 72-72. 

Kinks Kl and K2 perform annihilating random walks: 
Kl + Kl — > 0, K2 + K2 — ► 0, while the processes Kl + 
K2 — ► 0, K2 + K1 — ► are, however, forbidden. In other 
words, upon meeting, a Kl and a K2 "block" each other 
(do not annihilate and do not exchange sites) plj j. To 
our knowledge such kind of kinetics has not been studied 
before. Motivated by this fact we have decided to explore 
the critical behaviour of the above described system, on 
the line of compactness (q — 1), by computer simulation. 
In this study special attention will be paid to the p = 1/2 
symmetry point CDP2. 
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III. SIMULATIONS FROM RANDOM INITIAL 
STATE 

A. Kink decay simulations 

We have performed time dependent simulations start- 
ing from states with uniformly distributed species A, 
II, and 12, with respective densities: po(A), po(Il) and 
po(I2). At the CDP2 point unusual scaling behaviour 
of the density of kinks has been observed previously 
^lf : a deviation from the ordinary annihilation-diffusion 
process with kink-density decay p(t) ~ Instead, 

p(t) ~ t~ a with a « 0.55 has resulted from the first 
simulations. 

To check whether the observed deviation from stan- 
dard ARW behaviour is only a cross-over effect, or it 
heralds some basic feature of altered kinetics we have 
performed very long-time (t max = 10 6 MCS) simulations 
on systems with L = 24000 (Fig. |) Throughout 
the whole paper t is measured in units of Monte-Carlo 
sweeps. 



10 
time 



FIG. 3. Total kink number as the function of time 
started from symmetrical homogeneous random initial states: 
Po(/l) = Po(/2) = 0.3 (solid line) and p (Il) = p (I2) = 0.1 
(dotted line). The dashed line corresponds to the single 
species annihilating random walk (po(Il) — 0, po(I2) — 1/2) 
exhibiting p(t) oc t~ ' 5 . 

Figure |^ shows the results of simulations. It is seen 
that the deviation from the standard ARW value of the 
decay exponent remains present asymptotically as well: 
the local slopes of the decay curves 



a{t) 



In [p(t)/p(t/m)} 
ln(m) 



(1) 



(where we use m — 8 usually) go to constant values. 
Moreover, another interesting feature has become appar- 
ent: the kink-decay exponent depends on the initial con- 
centrations of the components po(Il) = Po(I2) and in 
such a way that for higher initial kink density (lower av- 
erage distance between the kinks) the decay is faster. 
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FIG. 4. Local slopes of the kink density decay for symmet- 
rical initial conditions po(H) = Po(/2) = 0.3,0.2,0.15,0.1,0.05 
(from bottom to top curves). The simulation result for one 
species (ARW) is also shown (top curve) 

Asymptotically, as po — ► 0, the average distance of dis- 
similar kinks goes to zero and the decay exponent tends 
to the ARW value: a — > 0.5. 

In the case of asymmetric initial condition (po(Il) =/= 
po(I2)) Kl-s and A2-s decay with different rates. The 
type that has smaller initial density decays faster. Ex- 
ample: in case of po(H) — 1/9, po{I2) = 1/3, K2 decays 
roughly like i~ a5 (unperturbed by Al-s) but the local 
slopes of the log-log p(Kl) — t dependence decrease from 
—0.5 strongly. 



IV. SIMULATIONS FROM AN ACTIVE SEED 

The cluster simulations |l^] were started from a state 
with uniformly distributed As and Jl-s except a single 
12 pair in the middle and the following characteristic 
quantities for the I2-s were followed: 

• the average number of /2-s, Ni2(t), 

• their survival probability P12 {t) , 

• and the average mean square distance of spreading 
of J2's from the center R 2 I2 {t) 

The above quantities were averaged over 7V S independent 
runs at the CDP2 point ( in case of R 2 I2 {t) only for sur- 
viving samples ). At the critical point we expect these 
quantities to behave for t — ■* 00, as 



N I2 (t) ex i" , 



P12H) cxr 



R 2 I2 (t) cx t 2 



(2) 



(3) 



(4) 
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Upon varying the initial density po(Il),ior the ex- 
ponents 5 and rj (defined similarly to eq.(||), the local 
slopes of Nj2(t) and Pniff) continuously changing val- 
ues have been observed (Figs. [|, [|). The deviation of 
these exponents from those of the single-species annihi- 
lation random walk process: 1/2 and 0, respectively is 
remarkable. The spreading exponent, z, on the other 
hand, seems to be constant within numerical accuracy 
and equals that of the single species annihilating random 
walk : z = 2/Z = 1 such that the generalised hyper- 
scaling law of the compact directed percolation p6| 



n + S = z/2 



(5) 



is satisfied. In this respect it is important that n has been 
found to be negative. 
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FIG. 5. Local slopes of the number of /2-s. The ini- 
tial state is uniformly distributed with initial densities: 
po(H) = 0.1 (solid line), 0.25 (dotted line), 0.5 (dashed line), 
0.75 (long-dashed line). 
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FIG. 6. The same as on Fig. |j for the cluster survival 
probability. 



V. CLUSTER SIMULATIONS OF COMPACT 
DIRECTED PERCOLATION CONFINED IN A 
PARABOLA. 

To understand the physics of our numerical results up 
to now we set up a parallelism with an other case where 
DP process is bounded by parabolic space-time bound- 
ary conditions. We perform simulations on the compact 
cluster version of this and compare the results with those 
of the GDK model in section VIII. 



Kaiser and Turban have investigated j23|,[24j the 1 + 1 
dimensional DP process limited by a special, parabolic 
boundary condition in space and time directions: 



V = ±ct k 



(0) 



where C changes under uniform length rescaling (by b) 
to: 



c =b Zk - 1 c 



(7) 



Here Z is the dynamical critical exponent. By referring to 
conformal mapping of the parabola to straight lines and 
showing it in the mean-field approximation Kaiser and 
Turban claim that for k < l/Z this surface gives relevant, 
for k > 1 /Z irrelevant and for k — 1 /Z marginal pertur- 
bation to the DP process. The marginal case results in 
C dependent non-universal power-law decay, (for details 
see next section), while for the relevant case stretched 
exponential functions have been obtained. The above 
authors have given support to this claim by numerical 
simulations. 

We have investigated the effect of parabolic and re- 
flecting boundary conditions for the CDP2 process nu- 
merically. Time-dependent cluster spreading simulations 
have been performed in the GDK model with parabolic 
boundaries such that at each time step the simulation 
region is bounded by two Jl-s at y m in and y m ax, where 



Vn 



Vn 



= L/2-2-C t k 



= L/2 



C t k 



(8) 



(9) 



Two 72-s have been put initially at the centre 
(L/2,L/2 + 1) and some initial space (two As to the 
left and right) between Kl-s and K2-s has been added. 
Therefore the role of Jl-s now is purely the formation 
of parabolic boundaries around /2-s and in fact we in- 
vestigate the plain CDP process with reflective boundary 
conditions. A typical 1+1 dimensional run looks like as 
shown below (1, 2 and stand for II, 12 and A, respec- 
tively): 
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When we fixed the exponent at k = 1/2, to make the 
situation marginal we found continuously changing ex- 
ponents for the exponents of the survival probability Pi% 
and the number of I2-s,Ni2, by varying the shape C (Fig. 

7) . One can see that the exponent-slopes of Ni2(t) (Fig. 

8) and those of Pi 2 (t) (Fig. |9|) change by varying C. The 
spreading exponent of the /2-s, z, seems to be constant: 
equal to unity. (Fig. These results are very similar 
to those of the seed simulations in GDK of section [V . 




FIG. 7. Parabola cluster confinement simulations for CDP. 
Middle curves: N I2 {t) (C = 2,1.5,1.2,1 top to bottom); 
Lower curves: P/2(t) (C = 2, 1.5, 1.2, 1 top to bottom); Upper 
curves: 

Rh(t) (G = 2 > 1-5, 1.2, 1 top to bottom). 
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FIG. 8. Parabola cluster confinement simulations for CDP. 
Local slopes of N12 for different values of C. Solid line: C = 2, 
dotted line: C = 1.5, dashed line: C = 1.2, long-dashed line: 
C = 1. 

The survival exponent changes in such a way that the 
hyper-scaling relation valid in case of compact directed 
percolation [|26|: 

z/2 = r] + 5 = 1/2 

is fulfilled. In this case it is important, again, that 77 
takes negative values. 
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FIG. 9. The same as Fig. g for the cluster survival proba- 
bility. 



The analysis based on local slopes (Figs, g, |ifj| ) 
shows again plateaus for high values of t, indicating true 
power-law behaviours. The magnitude of the exponent 
characterizing the decay of the density of J2-s decreases 
as C is increased reminiscent of a similar situation in [E3| . 
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FIG. 10. The same as Fig. [] for the exponent of the cluster 
spreading Rj 2 - 



VI. THEORETICAL CONSIDERATIONS FOR 
CDP CONFINED IN A PARABOLA 

A. Anisotropic scaling 

In 1+1 dimensional anisotropic systems the correlation 
length diverges as £y ~t~ u u in time and as £j_ ~t~" in the 
space direction with a dynamical exponent Z = vn jv (y is 
also denoted as v\_ in the literature). Covariance under 
a change of the length scales then requires two different 
scaling factors, b\\ = b z and b±=b . 

We consider now a system displaying anisotropic crit- 
ical behaviour and limited by a free surface in the (t, y) 
plane as given in eq(^). Under rescaling, with t' = t/b z 
and y' = y/b, C transforms according to eq.(0),as dis- 
cussed in the previous section. 

In the marginal case, which we will consider now, Z = 
1/k, the scaling dimension x m of the tip order parameter 
becomes C-dependent x m (C). 

The order parameter correlation function between the 
origin and a point at (t, y) transforms as 

G(A,t,y,i) =6- 2 ^G^A,^,|,^-) (10) 
when L is infinite. With b=t 1 / z , equation ( [To| ) leads to: 

Here lc = C 1 - Zk , A = ^ p ~ p ^ and x m is the scaling di- 
mension of the order parameter. The latter is connected 
to (3, the critical exponent of the order parameter via 
(3 = vx m . We will use this scaling form in the following. 
P is the usual order-parameter exponent, defined, for the 
DKCA, through p\ cx (p — p c )P , for p > p c \ pi is the 
stationary density of l's. In case of a first order transi- 
tion as is the case with compact directed percolation the 
following considerations hold. 



As already mentioned, 

P(t)oct~ s , (12) 

is the survival probability of l's for spread of particles 
(l's, in our notation) about the origin. Away from the 
critical point /?' governs the ultimate survival probability 
(starting from a localized source): P^ = lim^oo P(t) oc 

(p - Pcf ■ It is known that /3' = 1 in CDP. @ The order- 
parameter exponent, (3, however, is zero. This is because 
p = 1/2 marks a discontinuous transition, by symmetry. 
pi = for p < 1/2 and p\ = 1 forp > 1/2; strictly speak- 
ing, (3 is not defined here but it is natural to associate 
the value f3 = with the discontinuous transition. 

This problem with the ill-defined exponent (3 can be 
avoided following the lines of Grassberger and de la 
Torre's scaling argument |l2[ for discontinuous transi- 
tions as it has been presented by Dickman and Tretyakov 
. Consider a model with a transition from an absorb- 
ing to an active state at A = 0, with exponents 8, r\, 
z, and (3' defined as above. Suppose, however, that the 
order parameter p is discontinuous, being zero for A < 0, 
and 

p = p + /(A), (13) 

for A > 0, where po > 0, and / is continuous and van- 
ishes at A = 0. According to the scaling hypothesis for 
spreading from a source there exist two scaling functions, 
defined via Q 

p{y,t)~f- d *l 2 G{y 2 lt\^l^\ (14) 

and 

P{t) ~t" 5 $(At 1 / l 'ii). (15) 

(Here p(y, t) is the local order-parameter density, r ~ 
A _!y n.) Existence of the limit P^ implies that $(x) ~ 
x 13 as x — * oo, with (3' = Sun. In a surviving trial, the 
local density must approach the stationary density p as 
t — > oo, so p(y, t) ~ A^ po, for t — > oo with fixed y, and 
A small but positive. It follows that G(0, x) ~ x 13 for 
large x. 

An important consequence is that we can use as scaling 
dimension of the order parameter for CDP the value (3' 
in the relation x m — " instead of (3. Via scaling relations 
f3' = 6v\\, the values obtained by computer simulations 
for 8 will be compared with results for CDP+parabolic 
boundary conditions. In this context the connection, 
again via scaling relation, between 8 and the decay expo- 
nent of the density of kinks when starting from a random 
initial state a will also be made use of . 
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B. Mean field analysis for CDP confined in a 
parabola 

In this section we will follow the lines of the mean-field 
analysis of the 1 + I dimensional DP process confined 
by a parabola as given in ref. (^4|, but now applied to 
compact directed percolation. 

The basic processes are: 



On a parabolic system, we use the new variables t and 
C(t,y) =y/t k for which the free surface is shifted to £ = 
±C and equation is changed into 



dP 



1 d 2 P 



k C ^ 
dt 2t 2k d( 2 ' t d( 



(22) 



with the boundary condition P(t, ( = ±C) = 0. Through 
the change of function 



Qly+l 
Hy-i 



> > > 



P(t,0=exp[--C 2 t 2k - 1 }Q(t,C) 



P P 1 

FIG. 11. Site update rules for compact directed percola- 
tion. 

The order parameter correlation function is the prob- 
ability density P(t, y) for a site at (i, y) to be connected 
to the origin. 

First we consider the case without confinement. In 
mean-field approximation one can set up an equation for 
the connectedness at (t + l,y): 

P(t + l,y) = P {P(t, y+ 1)[1 - P(t, y-l)} + 

P(t,y-1) [ l-P(t,y+l)]} + P(t,y + l)P(t,y-l) (16) 

Going to the continuum limit the following differential 
equation is obtained: 



equation (E2|) leads to 



1 d 2 Q 



dQ = 

dt 2t 2k d( 2 



(k - l)?t 



t 



(23) 



Q (24) 



for which the variables separate when k = 1,1/2 or 0. 
These values of k just correspond to irrelevant, marginal 
and relevant perturbations. 

For k = 1 the critical behaviour is the same as for un- 
confined percolation as expected for an irrelevant pertur- 
bation. 

For the true parabola which is the marginal geometry, 
one may use equation (E2I) with fc = l/2 to obtain 



dP _ ld 2 P 



dt 2 dC, 2 



2~8C ^~t^ [ ' 



dP _ d 2 P 
dt ^ dy 2 



{2p-\)P+{l-2p)P 2 



(17) 



which is of the form studied in |g5[ for the directed walk 
problem. Writing Q(t,C) — (f>(t)ip(() leads to the follow- 
ing eigenvalue problem for i/KC) 



The homogeneous, stationary solution of eq.(|l7|) is: 



ld 2 ^ 
2~dC 



Cdip 2 
2dC 



(26) 



describing a first order transition for CDP at p c — 1/2, as 
it is the case. At the transition, p — p c , eq.(17) reduces 
to 



with (j)(t) ~t x . The solution is obtained as the eigen- 
value expansion 

oo 

P(t,y) = J2 B nt 



OP _ l d 2 P 
~dt ~ Ydy 2 



n=0 



x 2 1 - 

A «'2' 



2t 



(27) 



(19) 



This is the ordinary diffusion equation with RW solu- 
tion 



P(t,y) 



exp 



(-i) 



>2jvt 



(20) 



which is exact in the CDP case. From comparison with 
the scaling form in the previous subsection, the following 
(well-known) exponents for CDP arise: 



The behaviour at large t is governed by the first term 
in this expansion which decays as t~ °, i. e. with a C— 
dependent exponent as expected for a marginal perturba- 
tion. The dimension of the tip-to-bulk correlation func- 
tion is the sum of the tip and bulk order parameter dimen- 
sions, the first one being variable. Comparing with the 
form of the decay in eq.(|TT]) gives Ag = [x^/(C)+x m ]/Z 
and, using (21), the tip order parameter dimension is 
given by 



x™f(C) = 2X 2 



(28) 



v\\ = 1 v = 1/2 Z = 2 



_ 1 



(21) Its dependence on C is shown in figure 2. of 
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Analytical results can be obtained only in limiting 
cases which have already been discussed in [p5| . When 
C is infinite, A 2 , = 1/2, only the first term in the ex- 
pansion remains, which satisfies the initial and bound- 
ary conditions, giving back the free solution in equation 
(pp|) . For large C-values the tip exponent is x™f(C) — 

| + Cexp^— [1 + 0(e)] where e is the correction 

term itself. For narrow systems, the hyper-geometric 
function gives a cosine to leading order in C 2 . One ob- 
tains the following asymptotic behaviour in t 



the simple random walk and the random walk+exclusion 
(SEP) processes one also observes different dynamical be- 
haviours. This latter case is nothing else but the T = 
dynamics of the Id Ising Model with Kawasaki exchange, 
where we have different domain growth properties than 
in case of a simple random walk. 

An extensive numerical simulation with look-up table 
algorithm seems to confirm this expectation. As Figure 
|l2| shows the slopes of the density decay started from the 
special pairwise random states described above depend 
on the initial density p(0). 



and the tip exponent diverges as 7r 2 /4C 2 . 

For < k < 1/2 the dependence on t is expected to be 
a stretched exponential function. For details see p3|,p5[. 



VII. ANNIHILATING RANDOM WALK OF TWO 
SPECIES WITH EXCLUSION 



To check our results concerning the scaling properties 
of kinks in the GDK model at the CDP2 point we have 
carried out an explicit simulation of the annihilating ran- 
dom walk of two species (A, B) with exclusion. The 
model we have investigated has been suggested by Hin- 
richsen [^l] and is as follows. A(B) will hop to a neigh- 
bouring empty site with probability pi A (jalB) or annihi- 
late with a neighbour A (B) with probability p2A (p2B) 
while A and B do not react when getting into neigh- 
bouring positions. The initial configuration was chosen 
in such a way that allows pairs of the same kind to an- 
nihilate always within some finite time interval (i.e. the 
system evolves into an empty state), namely: 

... A.A.B.B.A.A.A.A.B.B 

That means that AA and BB pairs have been put in a 
Id ring with initial probability p(0). Had we not chosen 
the initial state like this the system would have ended up 
in some finite particle configuration where As and B- 
s follow each other alternatingly, separated by arbitrary 
empty regions. (This initial configuration is in agreement 
with the arrangement of the two kinds of kinks in some 
random initial state of the GDK model, too. ) The 
probabilities plA, plB, p2A and p2B have been chosen 
to be unity, to achieve maximum simulation effectiveness; 
no qualitative difference in the results have been found 
upon lowering them. 

Clearly this process is different from the simple an- 
nihilating random walk of two species A + B — ► Q, 
therefore we may expect that a field theory describing 
this model (which, however, is still missing) would result 
in a different fixed point with different critical exponents 
as well. Furthermore one can argue that when comparing 
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FIG. 12. Local slopes of particle decay in the annihilat- 
ing random walk + exclusion process of two species results 
(L = 24000). The initial conditions are: p = 0.1,0.25,0.5 
from bottom to top curve. 

The local slopes tend to constant values greater than 
a = 0.5 in agreement with the GDK kink results. The 
level-off in case of po = 0.5 happens only for t > 1.5 x 10 6 
MCS. The average AA and BB distances confining an 
other type of particle have also been measured during 
the simulations, that enables us to extract the amplitude 
(C) of the confinement in the function fitted (C x t~ a ). 
These values will be used to compare the results with 
those of the GDK (see next section) . 



VIII. SUMMARY OF TIME DEPENDENT 
RESULTS 

Since in all of the previously shown cases we found non- 
universal scaling depending on the initial conditions and 
the generic model to account for such behavior seems to 
be the CDP2 with parabolic boundary condition, we have 
decided to measure the region of confinement in all cases 
and plotted the survival probability exponents S and the 
kink decay exponents a as a function of the shape of the 
measured parabola. 

In the present case the final survival probability of a 
cluster plays the role of the order parameter exponent (3 
as explained at the end of Sec. VI and for the characteris- 



es 



tic exponents we have: S — (3 /i/±. Thus we have plotted 
the results for S(C). In a common graph the fitted values 
for a(C) are also shown; on the level of kinks the order 
parameter (3 is connected to a in the same way as (3' to 
S for 'spins' (see eg. p7[|). 

For random initial conditions in the GDK model the 
characteristic distance between two neighbouring kinks 
of a given type has also been measured. The average 
neighbour distance Ik2~K2 shown on Fig. [l^ have been 
obtained for initial densities: Pq{I1) = pa{12) = 1/3. 
The power law increase for large t (see the plateau for 
t > 30.000) with the same scaling exponent as the de- 
cay exponent is not very surprising because Pkink(t) cx 

1/IK2-K2- 

Since the K2 — K2 and the Kl — Kl pairs confine 
the motion of each other ( a Kl — K2 pair can not ex- 
change to K2 — Kl) this power-law increasing length 
scale imposes a 'stochastic' boundary condition (pres- 
sure on kinks) with a mean value of a parabola simi- 
larly that was investigated by Kaiser and Turban |f2~3f in 
case of 1 + 1 d DP processes |24| and adapted for the 
case of a CDP-like first order transition in section VI. As 
discussed before, the scaling dimension of the order pa- 
rameter changes continuously with the amplitude of the 
parabolically growing confining box size if it grows with 
the same exponent as the cluster inside. 



the results for GDK with random initial conditions: 




Po 


C 


a 


0.0 


oo 


0.5000(3) 


0.05 


14.09 


0.517(2) 


0.10 


7.62 


0.528(2) 


0.15 


6.11 


0.534(2) 


0.2 


5.85 


0.537(2) 


0.3 


4.18 


0.540(1) 



This is in agreement with Fig. 8 of (£4|, where an 
increasing C causes a decreasing exponent. Note that 
the first line in the table corresponds to the simple ARW 
process, therefore there is no confinement (amplitude C 
is oo). 

The main difference between the rigid parabola, bound- 
ary case and the "stochastic" conhnement is that in the 
latter case the boundary generates an additional noise to 
the motion of confined particles. Therefore we don 't sim- 
ple have "free" particles confined in a parabola, but they 
are also perturbed by the noise in such a Kl <-> K2 sym- 
metrical way that the outcome perturbation is marginal. 



Initial 11-12 distance 




ing size ~ boundary size= 



11-12 exert marginal 
perturbation on each other 



FIG. 14. Possible explanation for the non-universal scaling. 
The symmetric II an 12 clusters exert marginal exclusion 
perturbation on each other. 
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FIG. 13. Local slopes of the K2 — K2 neighbour distances 
in GDK model of size L — 24000. The initial state is uniform 
with po(il) = po(/2) = 1/3 

In our case we encounter similar situation. The kink 
density decay exponent a seems to vary continuously in 
case of symmetrical initial conditions. The initial condi- 
tions effect the amplitudes of the density decays (as was 
shown to be valid by field theory for pure the reaction 
diffusion of A,B particles [[|) and therefore the ampli- 



tudes of the confinement region sizes (C) (see Fig. 14). 
To compare our results with those of p3^4|] the form 
A + C t a has been fitted to the In2-K2\t) distances de- 
termined from the density decay simulations (assuming 
lR2-K2{t) — 2/ pkink(t)). The following table summarises 



In ARW simulations again the form A + C t a has been 
fitted for the measured AA, BB distances. 

In cluster simulations we fitted the form: y = C x t s for 
the region of confinements and determined the respective 
C-s in all cases. 
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-1.0 




FIG. 15. Confinement shape parameter (C) dependence of 
exponent 8 for cluster simulations in parabola (squares), 12 
cluster (diamonds) and a in simulations from uniform initial 
conditions in ARW2e (circles), GDK (triangles). The solid 
line shows the ir 2 /8C 2 mean-field approximation. 



As Figure 15 shows we obtained similar monotonically 
decreasing curves in all cases that also agrees with the 
results of ||§. The GDK-uniform and the ARW2 re- 
sults seem to lie on the same curve. The spreading sim- 
ulation results of GDK are different from those of the 
CDP+parabolic boundary condition case. 

This can be understood, however, since in the former 
case the confined particles (/2-s) have a back effect on 
the bulk (Jl-s) particles, while this is not the case when 
the boundary is fixed. 

We also show the 7r 2 /8C 2 curve, determined as the 
asymptotic solution C — > of the mean-field approxi- 
mation, see section VI eq.(^). This seems to be in fair 
agreement with the case of CDP+fixed parabolic bound- 
ary condition. 



IX. GENERALISATION FOR N > 2: 
SYMMETRIC ANNIHILATING EXCLUSION 
PROCESS OF N SPECIES 

We have carried out preliminary simulations in the 
generalized version of the model introduced in section 
VII. The system was started from configurations like: 

A..A.B...B..C.C..DD E . . EFF . . . 



where species of the same type can annihilate each other 
but different types can not exchange. Our results show 
that concerning the time dependence of the density the 
deviations from the square root decay persist for N > 2 
and this property seems to remain valid also for N — ► oo. 



20000 40000 60000 80000 100000 

time 



FIG. 16. Results for density decay in symmetric annihilat- 
ing exclusion process of N = 3, 4 (top, bottom curve) species. 



Fig. shows that a similar level-off can be observed 
in the local slopes as in the N = 2 case with a > 0.5. A 
tendency of increasing a with increasing N is apparent 
in our simulations; the growing effect of finite size cor- 
rections, however, prevented us from going further, for 
higher values of N, in this study. 



X. GDK ON THE LINE OF COMPACTNESS 

On the line q = 1 the role of the absorbing states (II 
and 12) is symmetric. The above reported simulation 
results for the GDK model refer to p = 1/2, the CDP2 
point. Now we will discuss the situation for p ^ 1/2. For 
p > 1/2 the creation of new As happens with probabil- 
ity greater than 0.5, the active domain size grows expo- 
nentially and the inactive regions die out quickly (and 
symmetrically); the all- A- phase plays the same role as 
the all- 0- phase of the original DK automaton 0. The 
deviation from the DK picture is quite apparent, however, 
for p < 1/2, as instead of the all-l-phase of DK, for all 
values of p with the exception of p — 0, Glauber-Ising - 
like kinetics governs the motion of kinks. The kinks here 
are extended objects (A-s), somewhat similarly to those 
in Grassberger's CA models ]iq] ( for p < p c , where p c is 
the critical point of a parity conserving phase transition 
) where also kinks of different extension (and there even 
of different structure) separate absorbing-phase clusters. 
(At p = diffusion stops and a striped space-time picture 
of II and 12 domains freezes in, again like in Grassberger's 
model A cited above.) The average size of As goes to 
zero between two domains of the same type quickly, the 
kink 'particles' of same type perform a biased random 
walk toward each other. On the other hand between do- 
mains of different types remains a film of As of average 
size 1, since a collision of an II and 12 domain always 
creates a new A at the next time step. That means that 
kinks of different types still block the motion of each 
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other. Therefore the role of As is similar to the kinks 
of the T = Glauber Ising model. On the whole line of 
1/2 > p > 0, in the long time limit, therefore one can 
expect the number of such kinks to decrease as (t~ 1 / 2 ). 
On the other hand this is a line of compactness, as all 
clusters growing from a seed are compact, the character- 
istic exponents (77, S and z) though strongly dependent 
on p and the composition of the initial state, satisfy the 
hyper-scaling law valid for compact cluster p6y . This 
statement involves that this line is a line of first order 
transition points with order parameter exponent (3 S = 0. 
This first order transition occurs in a symmetry-breaking 
('magnetic') field coupled to the II and 12 spins. For the 
detailed description of a similar situation see Jl8| ]. All 
these features have been supported by simulations. As 
an example we can give some results obtained at p = 0.4 
starting with a single 12 in the sea of II, As (25% of 
II, 75% of A ); S = 0.45, z/2 = 0.47, 7j J2 = -0.02. It 
is worth mentioning, that for hyper-scaling to hold it is 
again important that 77/2 is negative. 



XI. DISCUSSION 

We have investigated numerically the one dimensional 
generalized Domany-Kinzel cellular automaton on a line 
in the plane of its parameters where only compact clus- 
ters grow. The two types of kinks in the simplest ver- 
sion of this compact GDK model (two absorbing phases) 
follow annihilating random walk with exclusion (no reac- 
tion) between different types. The equivalence with an 
explicit two-species ARW model with exclusion is shown 
provided the initial state is prepared in such a way that 
the kinks are arranged in pairs with some density. High 
precision simulations revealed that this system relaxes 
in a non-trivial way: the decay exponent of the kink 
density depends on the initial density of kinks. We ar- 
gue that this is a kind of (internal) surface effect; sim- 
ilar to the ARW process confined by a rigid space-time 
parabola provided the power of the parabola is chosen 
to be marginal. This case has been explicitly investi- 
gated with the result that the spreading exponents be- 
have qualitatively the same way as expected from the cor- 
responding mean-field approximation. We have no proof 
of the marginality for the theory including fluctuations, 
but rely on symmetry arguments. If we assumed that 
particles would exert relevant perturbation (k < 1/2) 
on each other, the corresponding parabola picture would 
predict stretched exponential decay ( a behaviour that is 
very difficult to differentiate from power-laws by simula- 
tions) and the local slopes should go to some higher value 
as a function of time (meaning faster that any power-law 
decay). However our high precision data show just the 
opposite, the local slopes decrease as a function of time 
tending to a value somewhat greater than 1/2. 

Nevertheless, the possibility of pure square-root de- 
cay masked by some tremendously long crossover func- 



tion can not be ruled out. One could still expect a non- 
universal scaling of the survival probability of particles 
in the same way as was observed in 23 24j or in an other 
similar situation where a diffusing "prisoner" con- 
fined by marginally growing cage was investigated. In 
the latter case the the boundary condition was absorb- 
ing and an exact solution was possible giving an exponent 
for the survival probability which is a continuous function 
of the amplitude of the marginal parabola. Furthermore 
the survival of a diffusing prisoner (with diffusivity D) 
inside a cage where both walls diffuse (with diffusivity 
A) has been solved exactly and the decay exponent was 
found to be 7T /2 cos" 1 (D/(D + A) 

Non-standard scaling in a Id ARW model was also ob- 
served by Frachebourg et al. |30| . They have shown that 
the survival probability of particles in an ARW with one 
free boundary depends on the location of the particles. 
If we count the particles from the free boundary the sur- 
vival probability of odd particles decay with exponent 
0.225, while those of even numbered decay with expo- 
nent 0.865. The explanation for this is based on the fact 
that even numbered particles always have left and right 
neighbours during the process, while odd numbered par- 
ticles lack one of the neighbours and since the ARW in 
Id is diffusion limited they can escape. One can notice a 
similarity of this mechanism to the one in ARW2 models 
we investigated. Namely in our case there are infinitely 
many internal boundaries (generated by particles of dif- 
ferent types which can not exchange sites). 

Recently Bray |3l[] has shown that the relaxation to- 
wards the critical state in the 2d XY model depends 
on the initial state. This is very different from what is 
expected from field-theoretical RG predictions that can 
not take into account low-dimensional topological effects. 
Moreover, Bray has shown that the non- universal behav- 
ior of the persistence exponent in this case can be de- 
scribed by the random walk of a particle moving under an 
attractive central power-law force that creates marginal 
perturbation as compared to free random walk This 
scenario is similar to ours since we also have particles 
with RW exhibiting pressure of marginal strength on each 
other. 

Right after our submission two other preprints ap- 
peared on cond-mat p2| , ^3| , dealing with models very 
similar to those presented here and reporting results 
which are in accord with ours for those quantities they 
also investigated. 
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